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ABSTRACT 

A study of the blue compact dwarf (BCD) galaxy Mrk 996 based on high resolution 
optical VLT VIMOS integral field unit spectroscopy is presented. Mrk 996 displays 
multi-component line emission, with most line profiles consisting of a narrow, central 
Gaussian (FWHM^llO km s" 1 ) with an underlying broad component (FWHM^400 
km s _1 ). The broad Hi Balmer component splits into two separate broad components 
inside a 1".5 radius from the nucleus; these are attributed to a two-armed mini-spiral. 
This spiral-like nucleus rotates in the same sense as the extended narrow-line ionised 
gas but is offset by ^50 km s^ 1 from the systemic velocity of the galaxy. The rotation 
curve of Mrk 996 derived from the Ha narrow component yields a total mass of 
5x10 8 Mq within a radius of 3 kpc. From the Ha luminosity we infer a global star 
formation rate of ~2M Q yr _1 . 

The high excitation energy, high critical density [O in] A4363 and [N n] A5755 
lines are only detected from the inner region and exist purely in broad component 
form, implying unusual excitation conditions. Surface brightness, radial velocity, and 
FWHM maps for several emission components are presented. A separate physical anal- 
ysis of the broad and narrow emission line regions is undertaken. We derive an upper 
limit of 10,000 K for the electron temperature of the narrow line gas, together with 
an electron density of 170 cm -3 , typical of normal H n regions. For the broad line 
component, measured [O ill] and [Fe ill] diagnostic line ratios are consistent with a 
temperature of 11,000 K and an electron density of 10 7 cm -3 . The broad line emis- 
sion regions show N/H and N/O enrichment factors of ^20 relative to the narrow line 
regions, but no He/H, O/H, S/H, or Ar/H enrichment is inferred. Previous studies 
indicated that Mrk 996 showed anomalously high N/O ratios compared with BCDs 
of a similar metallicity. Our multi-component analysis yields a revised metallicity of 
>0.5 Z (12+log O/H = 8.37) for both the narrow and broad gas components, signifi- 
cantly higher than previous studies. As a result the narrow line region's N/O ratio is 
now typical for the galaxy's metallicity. The narrow line component's N/O ratio peaks 
outside the core region, spatially correlating with ~3 Myr-old stellar populations. The 
dominant line excitation mechanism is photoionisation by the ^3000 WR stars and 
~150,000 O-type stars estimated to be present in the core. This is indeed a peculiar 
BCD, with extremely dense zones of gas in the core, through which stellar outflows 
and possible shock fronts permeate contributing to the excitation of the broad line 
emission. 

Key words: galaxies: abundances - galaxies: individual (Markarian 996)- galaxies: 
kinematics and dynamics - galaxies: starburst - stars: Wolf-Rayet - galaxies: dwarf 
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1 INTRODUCTION 

Blue Compact Dwarf (BCD) galaxies provide a means of 
studying chemical evolution and star formation in low metal- 
licity environments in the nearby Universe. In general BCDs 
are faint (Mb > 18), have blue optical colours, are typ- 
ically small (<1 kpc) and are thought to be experienc- 
ing bursts of star formation in relatively chemically un- 
evo lved environments, ran ging from 1/2-1/50 solar metallic- 
ity l|Kunth fc Os tlin 2000]), and are therefore good analogues 
to high redshift star forming galaxies. 

The current study focuses on Ma rkarian 996, a BCD 
whose properties are far from ordinary. iThuan et al.l (|l996T l 
(TIL96 hereafter) presented a comprehensive HST study of 
Mrk 996 by analysing FOS UV and optical spectra and 
WFPC2 V- and /-band imaging. Their HST U-band im- 
age, shown in Figure [TJ revealed that the majority of the 
star formation in Mrk 996 occurs within a bright, compact 
nuclear H II region of ~3 arcsec angular radius (~315 pc 
for their adopted distance of 21.6 Mpc (TIL96)) surrounded 
by an extended elliptical low surface brightness component 
of ~40.6 arcsec diameter. However, the remarkable nature 
of the galaxy was exposed by the small-aperture (~1 arc- 
sec) FOS spectra obtained from the nucleus of Mrk 996. 
Line widths increase with the degree of ionisation; lines from 
singly ionised species being narrow and typical of H II re- 
gions whereas lines from ions such as 2+ and Ne 2+ are very 
broad, with FWHM widths up to ~900 km s" 1 (TIL96). Un- 
physically high [O ill] electron temperatures are derived for 
the nucleus unless an electron density (N e ) of at least ~ 10 6 
cm -3 is adopted, several orders of magnitude higher t han 
normal H II region electron densities (|lzotov et al.ll 19941) . In 
order to account for the observed line intensities, TIL96 ap- 
plied a two-zone density-bounded H II region model, with 
an inner core density of ~ 10 6 cm -3 and an outer zone with 
a density of ~ 450 cm -3 . They suggested that this large 
density gradient is caused by a mass outflow driven by the 
large population of Wolf-Rayet stars found to be present in 
the galaxy. 

Elemental abundances were also derived by TIL96, who 
estimated a metallicity of 0.22 Z®, based on the abundance 
of oxygen. As a result of their study, Mrk 996 was placed 
in a small group of eight BCDs k nown to display a sig - 
nificant overabundance of nitrogen (|Pustilnik et al.ll2004h . 
with derived N/O ratios that are 0.3-1.4 dex higher than 
for other BCDs of similar oxygen abundances. The en- 
hanced nitrogen with in this group has been suggested by 
iPustilnik et al. | (|2004h to be connected with merger events, 
in particular with a short phase of the consequent starburst, 
when many WR stars contribute to the enrichment of the 
interstellar medium (ISM). This was certainly the case for 
NGC 5253 (though t to be interacting with its companion 
galaxy NCG 5236 l|Moorwood fc Glass! Il982h ). where spa- 
tia lly resolved N / H and N/O enhancements were observed 
by I Walsh fc Rovl < |l987h to correlate strongly with WR star 
signatures. In support of this, a correlation was also found 
between spatial enhancements of He/H and N /O, attributed 
to he lium enrichment of the ISM by WR stars (|Walsh fc Rovl 
1 19891 ). 

The study of TIL96 showed that Mrk 996 has spatially 
varying physical properties. Spatially resolved kinematic and 
chemical abundance information across the galaxy is there- 



fore crucial for understanding the nature of this system; such 
information can be readily obtained with integral field spec- 
troscopy (IFS). In this paper we present high resolution op- 
tical IFS observations obtained with the VIMOS integral 
field unit (IFU) spectrograph on the 8.2m Very Large Tele- 
scope UT3/Melipal. The data afford us a new spatiokine- 
matic '3-D' view of Mrk 996. The spatial and spectral reso- 
lution achieved allows us to undertake a full multi-velocity- 
component analysis of this system, ultimately providing a 
more complete picture of its diverse ionised ISM. 

We adopt a distance of 22.3 Mpc for Mrk 996, at a 
redshifted velocity of 1642 km s _1 (corresponding to z = 
0.00544, thi s work), using a Hubble c onstant H = 73.5 km 
s _1 Mpc -1 jDeBernardis et al.ll200sh . 



2 VIMOS IFU OBSERVATIONS AND DATA 
REDUCTION 

2.1 Observations 

Two data sets were obtained with the Visible Multi-Object 
Spectrograph (VIMOS) IFU at the 8.2m VLT at ESO's 
Paranal Observatory in Chile. The VIMOS IFU consists of 
6400 (80x80) fibres coupled to microlenses and operates 
with four CCD detecto rs yielding four spect ral and spa- 
tial quadrant samplings (Zanichclli et al. 200^). The spatial 
sampling is contiguous, with the dead space between fibres 
below 10 per cent of the inter-fibre distance. The data sets 
were taken with the high resolution and high magnification 
(i.e. a partially masked IFU head containing only 40x40 fi- 
bres) settings of the IFU, resulting in a field-of-view (FoV) 
of 13" x 13" covered by 1600 spatial pixels (spaxels) at a 
spatial sampling of 0.33" per spaxel. Two different grisms 
were used, high resolution blue (HRblue, ~0. 51 A/pixel), 
covering 4150-6200 A and high resolution orange (HRor- 
ange, ~0. 60 A/pixel) covering 5250-7400 A; spectra from 
both grisms are illustrated in Fig [2] The instrumental width 
delivered by the IFU HRblue and HRorange grisms was mea- 
sured using arc lamp exposures. Gaussian profiles were fitted 
to a number of isolated lines within each extracted arc spec- 
trum and were found to have FWHMs of 2.3 ±0.1 A (113.5 
± 4.9 km s" 1 ) and 1.5 ±0.1 A (72.3 ± 4.8 km s" 1 ) for the 
HRblue and HRorange grisms, respectively. 

The observing log can be found in Table [T] Four expo- 
sures were taken per grism. The third exposure within each 
set was dithered by ±0.34" in RA and by ±0.58" in Dec 
(corresponding to a 2 spaxel offset in the X-direction) in 
order to remove any dead fibres when averaging exposures. 
All observations were taken at a position angle of ±30°. 

2.2 Data Reduction 

The data reduction was carried out using the GUI-based 
pipeline software GASGANcQ that allows the user to organise 
calibration files and run pipeline tasks. Three main tasks 
were used: (1) vmbias created the master bias frame; (2) 
vmifucalib determined the spectral extraction mask, wave- 
length calibration and the relative fibre transmission cor- 
rection; (3) vmifustandard created flux response curves from 

1 http:/ /www. eso.org/sci/data-processing/software/gasgano 
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Right Ascension 
Center: R.A. 01 27 35.54 Dec -06 19 35.1 



Figure 1. HST F569W WFPC2 image of Mrk 996 overlaid with the 13"xl3" VIMOS IFU aperture at a position angle of 30°, along 
with contours from the VIMOS IFU Ha narrow component emission line map (see Fig |3J. North is up and east is to the left. Ha 
contours are shown for the range 24.8 - 899.9, in steps of 119, in surface brightness units of xlO - 15 ergs cm -2 s — 1 arcsec - 2 . 



a summed spectrophotometric standard star spectrum. Fi- 
nally, the products of these tasks were fed into vmifuscience 
which extracted the bias-subtracted, wavelength- and flux- 
calibrated, and relative fibre transmission-corrected science 
spectra. Each VIMOS CCD quadrant is an independent 
spectrograph, therefore the IFU data processing was per- 
formed separately on each quadrant, creating four fully cal- 
ibrated 3-D arrays per science exposure. The flux calibra- 
tion was performed by multiplying the 2-D spectral frame 
for each quadrant by its relative response curve, derived 
from standard star observations in the same filter, grism 
and quadrant. The final data cube is produced by median- 
combining each exposure in a jitter sequence (i.e. re-aligning 
the dithered exposure with the other exposures of identical 
pointing) , after all the individual data quadrants for each ex- 
posure had been reduced. Reconstructing the cube spatially 
utilises the IFU table that lists the one-to-one correspon- 
dence between fibre positions on the IFU head and stacked 
spectra on the CCD (|Bastian et al.ll2006l ). A schematic rep- 
resentation of the dat a reduction processes can be found in 
IZanichelli et all l|2005l ). 

Unlike other IFU designs (such as VLT 
FLAMES/ Argus), there are no 'sky-dedicated' fibres 
in VIMOS IFU mode. Sky subtraction was performed by lo- 
cating a background region within each quadrant, summing 
the spectra over each of the spaxels in the reconstructed 
cube and subtracting the median sky spectrum of the 
region from its corresponding quadrant. Median combining 



was needed to ensure that any residual contamination from 
faint objects was removed. It should be noted that VIMOS 
IFU images are obtained as projections on the IFU CCD, 
i.e. they are mirror images to what is seen on the sky. 
Thus, all emission maps presented in this paper have been 
reflected about the Y-axis to obtain the correct orientation. 



2.3 Emission Line Profile Fitting 

The resolving power and signal-to-noise (S/N) of the data 
are high enough to allow the resolution of multiple compo- 
nents in the majority of the emission lines seen in the 1600 
spectra across the FoV. In order to fit spectra of this quan- 
tity and complexity, we utilised an automated fi tting proce- 
dure called PAN (Peak ANalysis; lDimedl2005h . This IDL- 
based, general-purpo s e curve -fitting package was adapted by 
IWestmoquette et all l|2007al ') for use with FITS format data 
by allowing multiple spectra to be read simultaneously in an 
array format. The user can interactively specify initial pa- 
rameters of a spectral line fit (continuum level, line peak 
flux, centroid and width) and allow PAN to sequentially 
process each spectrum, fitting Gaussian profiles accordingly. 
The output consists of the fit parameters for the continuum 
and each spectral line's profile and the x 2 value for the fit. 

It was found that in the spectra of Mrk 996 most high 
S/N ratio emission line profiles consist of a strong, narrow 
central Gaussian component (hereafter CI), superimposed 
on a very broad component containing a substantial pro- 
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Table 1. VIMOS IFU observing log 



Observation ID 


Date 


Grism 


Exp. time (s) 


Airmass range 


FWHM Seeing (arcsec) 


250520 


13/10/2006 


HR blue 


4 X 402 


1.135 - 1.097 


0.72 


250516 


19/10/2006 


HR orange 


4 X 372 


1.057 - 1.053 


1.15 




40 
U 35 

CM 

E 30- 



<" 25 
» 20 



15 



10 





IfwJ 



(b) Z 




5500 



6000 



6500 



7000 



Wavelength / A 



Figure 2. VIMOS IFU summed spectra of Mrk 996 after smoothing with a 5 pixel boxcar: (a) high resolution blue spectrum; (b) high 
resolution orange spectrum. Spectra are integrated over an area of 5.3x6.3 arcsec 2 (575x685 pc 2 ) and correspond to exposure times of 
4x402 s and 4x372 s for the HR blue and HR orange grisms, respectively. 



portion of the flux (hereafter C2). In addition to this, in 
the core of the galaxy the C2 component of the Balmer line 
profiles splits into two broad Gaussians of variable strength 
and velocity (hereafter C3 and C4, shifted towards the blue 
and red of the centroid of CI, respectively). The splitting 
of C2 into C3 and C4 and the region in which this occurs 
is discussed in Section [3.21 Fig [3] shows examples of fits to 
the Ha emission line from different regions of Mrk 996 in 
order to illustrate single, double and triple Gaussian fits, 
along with the corresponding fit residuals. In the same fig- 
ure a double Gaussian fit to the [O ill] A5007 profile is also 



shown corresponding to the same central spaxel for which 
Ha requires a triple Gaussian fit, thus illustrating the dif- 
ference between the Ha and [O ill] line profiles. Where ap- 
propriate, single or multiple Gaussians were fitted to the 
emission line profiles, restricting the minimum FWHM to 
be the instrumental width. Suitable wavelength limits were 
defined for each emission line and continuum level fit. Fur- 
ther constraints were applied when fitting the [S n] doublet: 
the wavelength difference between the two lines was taken 
to be equal to the redshifted laboratory value when fitting 
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the velocity component, and their FWHMs were set equal 
to one another. 

In order to rigorously determine the optimum number 
of Gaussians required to fit each observed profile the statis- 
tical F-test was used. The F-distribution function allows one 
to calculate the significance of a variance (x 2 ) increase that 
is associated with a given confidence level, for a given num- 
ber of degrees of freedom. We note that even though many 
low S/N ratio lines may in reality be composed of multiple 
emission components, it is often not statistically significant 
to fit anything more than a single Gaussian. This is why 
the robust F-test method was chosen for objectively deter- 
mining the optimum number of components to fit to each 
line. Resultant optimised fits for various line species are as 
follows: 

Single narrow Gaussian (i.e. CI component only) - 
[Fe in] AA4658, 4702, 4881, 4986, [O i] AA6300, 6364, [Nil] 
A6584, [S n] AA6717, 6731, and [On] AA7320, 7330. 

Single broad Gaussian (i.e. C2 component only) - [N n] 
A5755, [Om] A4363. 

Double narrow/broad Gaussian (CI and C2) - all 
Balmer lines (excluding the central spaxels), He I AA4471, 
5876, 6678, 7065, [O in] AA4959, 5007, [S in] A6312, and 
[Arm] A7136. 

Triple narrow/broad Gaussian fit (CI, C3, and C4) - all 
Balmer lines (central spaxels only). 

The errors reported by PAN during fitting underesti- 
mate the true uncertainties. We thus follow an error esti- 
mation procedure which involves the visual re-inspection of 
the line profile plus fit after knowing which solution was se- 
lected by our tests, and taking into account the S/N ratio 
of the spectrum. The errors on the fits are minimised across 
the IFU's aperture varying S/N ratio during the F-testing 
procedure by selecting the number of components which 
best-fit the profile within that particular spaxel. It is further 
noted that the uncertainties associated with each individual 
line component are coupled to those of the other compo- 
nents within a given line profile. By comparing PAN fits 
to those performed by other fitting techniques (e.g. iraf's 
SPLOT task) on line profiles with an established configura- 
tion (i.e. after performing the F-test) we find that uncer- 
tainties of ~5-10 and ~15-20 per cent are associated with 
the CI and C2-4 fits, respectively. A listing of measured flux 
and FWHM uncertainties appears in Table [5J where errors 
are quoted for individual component fits to emission lines 
detected on integrated spectra across the inner and outer 
regions of Mrk 996. 

2.4 Cube alignment and correction for differential 
atmospheric refraction 

Observations of an object's spectrum through the Earth's 
atmosphere are subject to refraction as a function of wave- 
length, known as differential atmospheric refraction, DAR. 
The direction of DAR is along the parallactic angle at 
which the observation is made. Each reduced data cube 
wa s corrected fo r this effect using the algorithm outlined 
by I Walsh fe Rovl (| 1990T ) : this procedure calculates fractional 
pixel shifts for each monochromatic slice of the cube rela- 
tive to a fiducial wavelength (i.e. a strong emission line), 
shifts each slice with respect to the orientation of the slit 
on the sky and the parallactic angle and re-combines the 



DAR-corrected data cube. Having data cubes with different 
wavelength ranges would ideally involve a choice of fidu- 
cial wavelength that is common to both, e.g. He I A5876 
for these VIMOS data. However, the S/N ratio of this line 
is not consistently high across the IFU aperture, and pixel 
shifts for each cube relative to this line were found to be un- 
reliable. Instead, we chose to correct each data cube using 
its strongest emission line: H/3 for the HRblue cube and Ha 
for the HRorange cube. The two cubes were then spatially 
matched by aligning the Ha and H/3 emission contours using 
a 2-D Gaussian fit. This involved a maximum shift of —0.73 
spaxels in X and +0.82 spaxels in Y. Two effects can intro- 
duce minor errors when aligning monochromatic H I images: 
(i) an extinction gradient across the observed region dispro- 
portionately affecting the two lines, and (ii) the presence of a 
temperature gradient, since the Hi lines have a very slightly 
different temperature dependence. However, both gradients 
would need to be significant for substantial errors to be in- 
troduced; the sub-spaxel shifts that were determined show 
that this is not the case here. 



3 MAPPING OF LINE FLUXES AND 
KINEMATICS 

3.1 Line fluxes and reddening correction 

Full HRblue and HRorange spectra are shown in Fig 
with identified emission lines labelled. Table|2]lists the mea- 
sured FWHM and observed and de-reddened fluxes for the 
main velocity components of the detected emission lines. 
The fluxes are from spectra summed over the core region 
and the outer annulus region, as defined in Section 13.41 
and are quoted relative to the flux of the corresponding 
H/3 component. They were corrected for re ddening using the 
Galactic reddening law of iHowarthl l|l983h using c(H/3) val- 
ues derived from the Ha/H/3 and H7/H/3 line ratios of their 
corresponding components, weighted in a 3:1 ratio, respec- 
tively, in conjunction with the theoretical Case B ratios from 
iHummer fe Storevl (|l987T ). Following the same method, an 
average c(H/3) map was also derived using ratioed Ha/H/3 
and H7/H/3 emission line maps. The c(H/3) map was used 
for correcting emission line maps for reddening in order to 
then create abundance maps. In the direction to Mrk 996 a 
foreground Milky Way reddening of E(B — V) = 0.04, cor- 
respond ing to c(H/3) = 0.06, is indicated by the extinction 
maps of ISchlegel et all (|l998l ). Values of e(H/3) of 0.10±0.02 
and 0.51±0.08 are applicable to the narrow component emis- 
sion in the inner and outer regions of Mrk 996, respectively, 
while values of 0.42±0.04 and 0.91±0.12 are applicable to 
the broad component emission for the inner and outer re- 
gions, respectively. 

3.2 The Balmer line properties 

All the Balmer lines detected from Mrk 996 (Ha, H/3 and 
H7) show a common multi-component velocity structure 
which in the innermost spaxels is exclusive to the Balmer 
lines, i.e. it is not seen in the forbidden lines. Here we present 
a discussion of the highest S/N ratio Hi line, Ha. Maps of 
the individual Ha velocity components (flux, radial velocity 
and FWHM) are shown in Fig [4] As mentioned previously, 
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Figure 3. Example Mrk 996 emission line profiles and Gaussian fits: (a) a typical spaxel from the core region with a triple Gaussian fit 
to the Ho line; (b) the same spaxel spectrum shown in (a) but with a double Gaussian fit [note the larger residuals of the fit as compared 
to (a)]; (c) Ho emission line from a spaxel in the outer region which required a double Gaussian fit; (d) Ho emission line from the edge 
of the outer region where a single component fit is sufficient; (e) the H/3 emission line for the same core spaxel as that in (a), displaying 
the same triple component structure (H7 also shows an identical structure in the core spaxels); (f) the [Oin] A5007 emission line for the 
same core spaxel as that of (a) and (b), illustrating an optimal double Gaussian fit (rather than the triple components required for the 
Balmer lines). The flux units are xl0 — 15 erg s — 1 cm" 2 A -1 . The fit residuals are plotted under each panel (in units of cr) with dashed 
guidelines at ±1<t. The wavelengths are not corrected for redshift. 
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Figure 4. Maps of Mrk 996 in the Ha velocity components CI— 4; (a) logarithmic surface brightness; (b) radial velocity (relative to the 
heliocentric systemic velocity of +1642 km s — 1 ), where the white dotted line overlaid on the CI velocity map represents a cut along the 
rotational axis, as used for the position- velocity diagrams shown in Fig [5] (c) FWHM corrected for the instrumental PSF. See text for 
details. 
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Table 2. Emission line measurements for summed spectra of the inner core (1.7"x2.3") and an 'outer annulus' region whose extent is 
5.3"x6.3" but from which the inner core spectrum has been subtracted (see Section l3,4l l, Observed fluxes and dereddened intensities are 
given for the main two velocity components and are relative to H/3 = 100 (the corresponding velocity component of H/3). The employed 
logarithmic extinction constants, c(H/3), are quoted (see Section l3.1l l; integrated H/3 fluxes for the two main velocity components, in units 
of xlO 15 erg cm -2 s _1 , are also listed. 
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6300 [Ol] 


CI 


128. Ozt 13.8 


4.6± 0.4 


4.3 zt 0.4 


141. 9zt 22.9 


7.6zb 1.5 


5.4 zt 0.7 


6312 [Sm] 


CI 


107.4± 13.8 


1.2± 0.2 


1.2 ± 0.2 


133.0zt 31.8 


3.8zh 0.8 


2.7 zt 0.4 


6312 [Sm] 


C2 


403.0zt 29.0 


5.8± 0.4 


4.4 ± 0.3 


_ 


_ 


_ 


6363 [Oi] 


CI 


153.2± 20.3 


1.7± 0.2 


1.6 zt 0.2 


147.5zt 38.2 


4.0zb 1.1 


2.9 zt 0.6 


6563 Ha 


CI 


116.5zt 0.9 


306.9± 4.4 


285.6 zt 3.6 


118. 8zt 0.4 


418.9zh 65.3 


287.6 zt 1.6 


6563 Ho 


C2 


420.5zt 7.3 


400.5± 7.5 


294.0 zt 4.6 


407.7zt 9.6 


555. 5zb 85.3 


284.2 zt 0.8 


6584 [No] 


CI 


126.2± 4.1 


38. 0± 1.1 


35.3 zt 1.0 


123. 4zt 2.2 


45.4zb 7.1 


31.1 zt 0.4 


6678 He I 


CI 


112.7zt 4.9 


3.5± 0.2 


3.2 zt 0.2 


124.8dz 14.8 


6.8zb 1.3 


4.6 zt 0.5 


6678 He I 


C2 


550.3zt 28.8 


9.2± 0.4 


6.6 zt 0.3 








6717 [Sn] 


CI 


118.8± 1.8 


30.9± 0.5 


28.6 zt 0.5 


117.4zt 2.2 


31.6zb 5.0 


21.2 zt 0.4 


6731 [Sn] 


CI 


118.5zt 1.8 


24.6± 0.5 


22.8 zt 0.4 


117.2zt 2.2 


22.4zh 3.5 


14.9 zt 0.4 


7065 He I 


CI 


76.00zt 65.4 


1.1± 0.9 


1.0 zt 0.8 


217.8zt 56.1 


7.0zb 1.9 


4.4 zt 1.0 


7065 He I 


C2 


425.4± 58.2 


24.4± 2.8 


16.8 zt 1.9 








7136 [Arm] 


CI 


127.3zt 17.2 


10.3± 2.3 


9.4 zt 2.1 


127.8zt 11.4 


13.5zh 2.3 


8.5 zt 0.7 


7136 [Arm] 


C2 


441. Orb 213.9 


8.2± 3.5 


5.6 zt 2.4 








7320 [On] 


CI 


176.6zt 11.9 


6.1± 0.3 


5.6 zt 0.3 


219. 2zt 62.3 


9.2zb 2.3 


5.6 zt 1.1 


7330 [On] 


CI 


165.3zt 11.9 


5.8± 0.3 


5.3 zt 0.3 


111.3zt 21.3 


10.8zh 2.2 


6.6 zt 0.8 


c(H/3) 


CI 




0.10±0.02 






0.51zh0.08 




c(H/?) 


C2 




0.42±0.04 






0.91zh0.12 




F(K0) 


CI 




626.30±4.31 






180.33zh28.13 




F(K{3) 


C2 




431.70±4.46 






36.67zh5.57 





up to three components are resolved in the Ha line profile 
in the nuclear region of Mrk 996: 

(i) A central, narrow Gaussian (CI) is detected through- 
out the galaxy (see the flux map in Fig peaking at 
a = 01 h 27 m 35!5±0.1, 5 = -06° 19' 36'.'2±0.2 (J2000), with 
a surface brightness of 10.1 xl(F 12 erg cm~ 2 s _1 arcsec -2 
in the central region, and gradually decreasing over a region 
with a major axis radius of 3". 84 (403 pc) and a minor axis 
of 2". 39 (251 pc); beyond this the Ha emission traces a very 
low surface brightness area (which TIL96 fitted with an ex- 
ponential disk on their HST images). The galaxy does not 
show a smooth elliptical shape when viewed in Ha. Fig [T] 
shows the flux contours of the narrow Ha component (CI) 
overlaid on the HST F569W ( F-band) image. The Ha emit- 
ting region imaged by VIMOS is mostly confined within the 



inner part of the F569W image but the major axes of the 
two are not aligned. A dust extinction feature is present in 
the north-west section of the inner star forming region that 
distorts the F569W flux isophotes (TIL96). The dust fea- 
ture also distorts the outermost Ha contours of the VIMOS 
data. 

(ii) In the very core of Mrk 996, over an area of ~l".5 
in radius (162 pc), the narrow CI profiles are accompanied 
by broad underlying components that are best fit with two 
Gaussians (C3 and C4) to the blue and red of CI, respec- 
tively (Fig 3^,). The significance of these features is discussed 
below. 

(iii) Beyond a radius of ~l".5 from the core, C3 and 
C4 blend into a single broad Gaussian (C2) (Fig |4ji). The 
region over which this component is seen is ~2-4 arcsec in 
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Figure 5. Position velocity diagrams for He*: (a) along the axis of rotation, defined as 5° east of north, as indicated by component CI 
in Fig |4p; (b) along the direction 90° to the rotation axis. Both diagrams are relative to a systemic velocity of +1642 km s — 1 . Typical 
radial velocity errors are <3, 8—20, 12—15, and 7—11 per cent for CI, C2, C3 and C4, respectively, with a spatial error of 0.5 spaxels 
(±0.17"). 



radius (216-432 pc); its outer boundary matches well that 
of the brightest portion of the CI flux map. 

The radial velocity map of component CI (Fig 2p) 
shows a strong velocity gradient from west to east, with 
a definable axis of rotation at PA = 5° (marked by the 
dashed line across the CI velocity map). Neither the rota- 
tion axis as revealed by VIMOS nor the axis normal to it 
are aligned with the apparent, i.e. photometric, major axis 
of the galaxy which lies at PA ~ 23° (TIL96). A velocity gra- 
dient in approximately the same direction is observed in the 
broad component C2, although in this case the negative ve- 
locity component is not exactly balanced by a positive coun- 
terpart. Fig [5]shows position- velocity (P-V) diagrams along 
the adopted axis of rotation (Fig [5ji) and also perpendicular 
to this axis (Fig Od). A new estimation of the heliocentric 
systemic velocity of Mrk 996 can be made using Fig [5^,; a 
radial velocity of 1642 ± 10 km s _1 is needed to normalize 
the distribution of CI velocities along the axis of rotation to 
zero (cf. the 1622 ± 10 km s _1 estimate of TIL96). The nar- 
row line component CI shows a velocity gradient along the 
direction normal to the rotation axis, with radial velocities 
ranging from ~ —50 km s^ 1 to +50 km s -1 , and the P-V 
diagram of component CI shows a linear velocity gradient 
indicating solid-body rotation (Fig [5jj) . This type of struc- 
ture is not uncommon for BCDs observ e d in t he Hi 21 cm 
line and has been used by Ivan Zee et al.1 (|200ll ) to constrain 
evolutionary scenarios. 

Assuming a Keplerian velocity profile we have converted 
each radial velocity data point of component CI (Fig [5p) 
into a mass within radius r, as shown in Fig HJ] This plot 
shows dynamical mass growing at similar rates on either 
side of the centre, rising to a total mass of ~ 5±lxlO 8 M0 
within a radius of ~0.75 kpc. This is in good agreement with 
the total mass of 4.3 x 10 8 M q estimated by TIL9 6 from H I 
emission features detected by Thuan et al.l fl995). It should 



be noted that all previous kinematical s tudies of BCDs have 
been based on Hi 21cm observations l|van Zee et al.ll200ll ; 
iThuan et aT1l2004r i. 

The broad line component C2 also shows a velocity gra- 
dient roughly perpendicular to the rotation axis, (Fig 
aligned more closely with the major axis of Mrk 996, with 
radial velocities that merge with those of CI 2" from the 
centre of the galaxy, but extending to larger negative val- 
ues. No definite velocity structure can be seen along this 
axis for the C3 and C4 components; however, these show an 
opposing velocity gradient along the galaxy's rotation axis 
(Fig [SJi). Their radial velocities lie much higher at —60 to 
-160 km s _1 for C3 and +110 to +210 km s" 1 for C4, with 
the gradients in both showing a similar 'S' shape offset by 
1" from the dynamical centre of Mrk 996. It is further noted 
that the velocity centre of symmetry of the 'S' feature is 
offset by +50 km s _1 from the heliocentric systemic veloc- 
ity of Mrk 996 measured above. The 'S' kinematic feature 
can be attributed to a two-arm spiral structure located at 
the nucleus of Mrk 996, whose approaching and receding 
arms are respectively traced by the velocity distribution of 
the C3 and C4 components of Ha. Furthermore, our anal- 
ysis reveals that the projected angular velocity vectors on 
the plane of the sky of the spiral and of the gas traced by 
the CI narrow line component are at an angle of <90°. In 
their WFPC2 images of Mrk 996, TIL96 identified a small 
spiral structure along the east-west direction (see their fig 
4b) whose pivot is the nuclear star-forming region which 
is slightly offset from the centre of the outermost F791W 
image isophotes. The pivot of the spiral and the peak flux 
position on the VIMOS Ha CI map coincide spatially. The 
1" offset mentioned above for the spatial centre of symmetry 
of the 'S' feature (Fig [5^,) is consistent with the asymmetry 
noted by TIL96. This fact, along with the ~+50 km s _1 
offset of the spiral's pivot with respect to the systemic ve- 



10 B. L. James et al. 



Distance from FoV Centre / arcsec 

2 4 6 



0.0 0.2 0.4 0.6 

Radius / kpc 

Figure 6. Dynamical mass as a function of radius, r, derived 
from the distribution of CI radial velocities along the direction 
normal to Mrk 996's rotation axis: distance east from the centre 
represented by crosses and distance west from the centre repre- 
sented by stars. 



locity of the galaxy, indicates that the spiral-like nucleus of 
Mrk 996 is kinematically decoupled to some degree from the 
main ionised gas component. This could indeed be the fossil 
kinematic signature of a past merger event. 

The FWHMs of Cl-4 in Ha are shown in Fig \4jp, cor- 
rected for the instrumental PSF (see Section \2. 2 \ . A peak in 
the FWHM of the narrow CI component, ~110 km s _1 , is 
seen at the nucleus of Mrk 996. This gradually, but not uni- 
formly, decreases to ~65 km s _1 towards the outer regions. 
The gradient in the FWHM of C2 is far stronger, peaking 
at ~450 km s _1 at two positions northwest and southest 
of the nucleus and decreasing to ~100 km s _1 elsewhere: 
the broader line widths along the former direction could be 
the signature of an expanding shell surrounding the inner 
star forming region, blown outwards by the stellar winds of 
the core cluster. Higher resolving power IFU spectroscopy 
is needed to disentangle the kinematics of this region. No 
significant gradient is seen in the FWHM of components C3 
and C4, both remaining between 200-300 km s _1 . 



3.3 The [O ill] nebular line properties 

As mentioned previously, the structure of the H I Balmer 
lines is not mirrored in the forbidden nebular lines. All the 
strong forbidden lines (with the exception of [O ill] A4363 
and [N n] A5755 discussed below) are comprised of a nar- 
row, central Gaussian (CI) with an underlying single broad 
Gaussian component (C2). An optimal fit, after applying 
the rigorous statistical F-test to a spaxel from the core of 
Mrk 996, can be seen in Fig |3}'. For the spaxels in which the 
Balmer components C3 and C4 are detected, the FWHM of 
the [O ill] A5007 C2 component is ~390-415 km s" 1 and 
that of each of Ha C3 and C4 is ~180-240 km s _1 . The oc- 
currence of additional C3-4 components in the Balmer lines 
is in contrast to the non-occurrence of C2 splitting in the 
[O ill] nebular lines and can perhaps be ascribed to a lower 
S/N ratio for the broad component of [O ill] AA4959, 5007 
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Figure 8. Flux distribution of [O III] A4363 per 0.33x0.33 arcsec 2 
spaxel with Ha contours overlaid. The area over which A4363 is 
detected is used to define an 'inner core' region, with a surround- 
ing 'outer region' covering the full extent of the Ha emission. 



compared to the C2 of H I lines: the observed C2/C1 flux ra- 
tios pertaining to the [O ill] nebular lines in the inner core of 
Mrk 996 (as defined in the following subsection) are half of 
those observed in the Hi lines (Table 2). However, from our 
analysis of the broad C2 component, we find log (iV e /cm~ 3 ) 
= 7.25 ± 0.25 from the [O ill] lines (Section gj. This is high 
enough to collisionally suppress AA5007+4959 whose critical 
density for de-excitation is 6.4xl0 5 cm -3 . The distribution 
of flux in the CI and C2 components of the [O in] A5007 
line is shown in Fig [7] The shape and extent of the emitting 
regions are similar to those of the corresponding Ha CI and 
C2 components. 



3.4 Broad emission lines: [O in] A4363 and [Nil] 
A5755 

Our multi-component line analysis shows that only the 
broad (C2) components of the auroral [O ill] A4363 and [N n] 
A5755 transitions are present, and only in the inner galaxy. 
The implications of this for the use of [O ill] A4363 as a 
temperature diagnostic are discussed in Section [4] For the 
purpose of this study the A4363 spectral map (Fig [8} is used 
to define a two-region model of Mrk 996: an 'inner core' re- 
gion, defined by A4363 emission extending over an area of 
1.7x2.3 arcsec 2 (185x250 pc 2 ), and an 'outer region' large 
enough to incorporate almost the full extent of Ha emission 
over an area of 5.3x6.3 arcsec 2 (575x685 pc 2 ). The lack of 
any detectable [O ill] A4363 emission in the outer region is 
illustrated in Fig [5] Here we compare the summed spectrum 
of the outer region plus core with the summed spectrum of 
the core region, and also show the residual outer region spec- 
trum. An emission line map of the broad [N n] A5755 line 
is shown in Fig 110b .; its flux distribution correlates strongly 
with that of the similarly broad A4363, suggesting that the 
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Figure 7. [Om] A5007 emission line maps of Mrk 996 showing the flux distribution per 0.33x0.33 arcsec 2 spaxel: (a) in the narrow CI 
component; (b) in the broad C2 component. Overlaid are the Ha narrow component contours (white), as shown in Fig [T] 



700 



150 









- Sum lnner+ Outer Region _ 




Sum Inner Region 




Outer Region _ 








I [0III]X4363 





4300 4320 4340 4360 4380 4400 

Wovelength (A) 

Figure 9. Summed spectra from the inner core and the outer 
region illustrating the lack of detectable [O III] A4363 in the outer 
region of Mrk 996. 



excitation of these lines is affected by processes that are par- 
ticular to the inner core region of Mrk 996. In Fig [11] the 
velocity profiles of the broad C2 components of Ha, H/3 and 
[O ill] A5007 are shown, overlaid in velocity space with the 
profiles of the [Om] A4363 and [Nil] A5755. Each spectral 
line has been summed over the inner core region and all 
show a FWHM of ~300 kms -1 , suggesting that a common 
excitation mechanism is responsible for the emission of the 
broad component. The larger widths of [O ill] A4363 and 
[N n] A5755 indicate that the physical conditions within the 
inner core are different from those in the outer region. 
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Figure 11. Velocity profiles of the broad C2 components of Ha, 
H/3, [O m] A5007, and of [O in] A4363, [N n] A5755 emission lines 
(relative to the heliocentric systemic velocity of +1642 km s — 1 ). 
Spectra have been summed over the inner core 1.7x2.3 arcsec 2 
region of Mrk 996. 



4 MULTI-COMPONENT DIAGNOSES: 
ELECTRON TEMPERATURE AND 
DENSITY 

The multi-component nature of the emission lines from 
Mrk 996 must be taken into account when applying plasma 
diagnostic methods. For example, since [O ill] A4363 is only 
detected from the central region of Mrk 996 and completely 
lacks the narrow line component that [O m] AA4959, 5007 ex- 
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Figure 10. Flux distribution per 0.33x0.33 arcsec 2 spaxel of (a) [Nil] A5755 for which only a broad component is detected and which 
is seen only within the core region; (b) [Nil] A6584, whose flux distribution peaks in the core and towards the north-east; a fainter peak 
is seen west of the core. Overlaid are the Ha narrow component contours as shown in Fig [T] 



hibit (see Fig EJ, an integrated [O ill] (A5007+A4959) /A4363 
ratio can not provide a useful temperature diagnostic; this 
would lead to incorrect results. Previous analyses lacked the 
spectral and spatial resolution needed to decompose the line 
profiles and were forced to assume a single density of N e 
~10 6 cm -3 in order to match the anomalously low inte- 
grated [Om] (A5007+A4959)/A4363 flux ratio (TIL96). 

The detection of separate narrow and broad line compo- 
nents opens up the possibility of determining separate elec- 
tron temperatures and densities and conducting a separate 
abundance analysis for each. 



4.1 The Narrow Component Gas 

In order to estimate the electron temperature (T e ) from 
which the narrow line component CI emission arises the fol- 
lowing method was used. An upper limit to the T e applicable 
to the H I and [O ill] nebular CI components can be obtained 
via the simulation of a narrow [O ill] A4363 component using 
the observed narrow A5007 component. Predicted A4363 CI 
line intensities were obtained by multiplying the observed 
A5007 CI profile by the theoretical A4363/A5007 intensity 
ratio for a range of T c 's at iV c = 170 cm -3 ; the latter was 
measured from the [S n] doublet ratio as described below. It 
should be noted that the [S n] lines only show narrow CI 
components (Table . The scaling factors were applied to 
both the inner core and outer region CI component A5007 
absolute fluxes. Simulated narrow A4363 profiles of decreas- 
ing intensity were then subtracted from the observed A4363 
line until the residual was no longer detectable, as shown 
for the core region A4363 profile in Fig 1121 It is found that 
a measurable narrow A4363 CI component would be de- 
tectable at an electron temperature of > 10,000 K for both 
the inner core and the outer region. Thus we adopt an up- 
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Figure 12. Inner core region dereddened [Our] A4363 profile 
(solid red curve) with simulated narrow components subtracted 
from it, corresponding to the dereddened [Om] A5007 narrow 
component profile, normalised to the different A4363 flux levels 
that correspond to different adopted T e values. The subtracted 
narrow component profiles correspond to electron temperatures 
of 9000 K, 11,000 K and 13,000 K for an electron density of 170 
cm -3 . An upper limit of 10,000 K is estimated for the T c of the 
gas emitting the narrow CI component in the inner core region. 
See text for details. 



per limit of 10,000 K to the electron temperature of the gas 
emitting the CI component throughout Mrk 996. In Sec- 
tion [52] we argue that this limit is closer to 9200 K, based 
on the inference that these O/H abundance ratio across the 
narrow and broad line emitting regions is invariant. 

The [S n] AA6717, 6731 lines, whose intensity ratio is a 
common electron density diagnostic for H II regions, show 
only the narrow emission component and were used to com- 
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Table 3. Electron densities derived from the narrow line [Sn] 
A6716/A6731 intensity ratios (1.26±0.03 and 1.41±0.03, respec- 
tively for the inner core and outer region of Mrk 996). 



Adopted T e (K) 


N c (cm" 3 ) 


N c (cm" 3 ) 




Core Region 


Outer Region 


9,000 


170±40 


10± 30 


11,000 


170±40 


10± 30 


13,000 


170±40 


10± 30 



pute density values for a range of electron temperatures for 
the inner and outer regions of the galaxy (Table [3] using the 
iRAF'fl temden task of the nebula package). These densi- 
ties are representative of conditions in the gas from which 
the narrow component emission arises. This propels us to 
use alternative methods to investigate the density of the gas 
emitting the broad line components in the inner core region 
of Mrk 996. With the electron temperature of T c < 10,000 K 
obtained above, we find that electron densities of 170 cm -3 
and 10 cm -3 are representative of the narrow CI component 
gas in the inner and outer regions of Mrk 996, respectively. 
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Figure 13. Theoretical [Fein] line flux ratios, for A4881 (red) 
and A5270 (black), relative to A4658 as a function of electron den- 
sity for temperatures between 7,000 - 20,000 K. The horizontal 
lines denote the observed flux ratios (solid) and their uncertain- 
ties (dashed) from the summed spectra of the inner core region 
of Mrk 996. 



4.2 [Fein] line diagnostics 

[Fe in] emission line ratios provide useful density diagnos- 
tics th^_cari_trace_both high and low electron densities; 
see iKeenan et al.l ((2001). The summed inner region spec- 
trum of Mrk 996 (Fig [5]) shows a few of the brightest 
transitions among the 3d 6 levels of Fe III; [Fein] AA4658, 
4702, 4881, 4986, and 5270. The flux in these lines above 
the continuum level was measured by integrating over the 
line profile. To maximise accuracy, the fluxes for A4658 
and A4702 were measured after removing any contamina- 
tion from Wolf-Rayet stellar features by subtracting the WR 
template fit derived in Section [6] Fig [13] shows the extinc- 
tion corrected fluxes of A4881 and A5270 relative to A4658 
along with theoretical ratios a s a function of N e , for T c 's 
between 7,000-20,000 K, from IKeenan et al.l (|200lh . Inter- 
polating between the theoretical data for an adopted T c of 
10,000 K, the A4881/A4658 ratio indicates log (N c /cm~ 3 ) = 
5.7 ± 0.3, which falls within the somewhat larger range of log 
(iV c /cm- 3 ) = 6.5±°{ B 5 implied by the observed A5270/A4658 
ratio. These results indicate the presence of a very dense 
ionised medium throughout the inner core of Mrk 996. 

In contrast to the above however, the [Fein] 
A4986/A4658 ratio is a more sensitive tracer of low electron 
densities. Fig [l4]shows that relative to A4658 the A4986 line 
should not be detectable at densities higher than ~ 1,000 
cm -3 ; the fact that the line is detected from Mrk 996 indi- 
cates the presence of a low density gas component as well. 
This is not inconsistent with the high densities derived above 
for the inner core as the spatial distribution of A4986 emis- 
sion (not shown) contrasts strongly with that of the other 
[Fe ill] lines - rather than peaking in the core of the galaxy, 
A4986 is distributed throughout Mrk 996 and is not con- 
fined to the inner region defined by the broad [O ill] A4363 

2 I RAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy 
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Figure 14. Theoretical [Fein] line flux ratios, for A4986 (black) 
and A4207 (red), relative to A4658 (in the region in which A4986 is 
seen) , relative to A4658 as a function of electron density for tem- 
peratures between 7,000—20,000 K. The horizontal lines denote 
the observed flux ratios (solid) and their uncertainties (dashed) 
from the summed spectra of the region where A4986 is detected. 



emission. There is a clear detection of A4986 emission over a 
group of 8 spaxels extending ~1".6 east from the core. The 
A4986/A4658 and A4702/A4658 intensity ratios integrated 
over these 8 spaxels (Fig I14[) confirm that this is a low den- 
sity area (iV G ~100-1000 cm' 3 ), typical of normal H II re- 
gions. This emission could be tracing lower density ionised 
gas throughout Mrk 996. In a st udy of high excitat ion neb- 
ulae in the Magellanic Clouds bv lNaze et al.l l|2003h . [Fe III] 
A4986 was only detected in the shell of the super-bubble 
surrounding the bright Ha nebula N44C. Due to the large 
distance of Mrk 996, we are unable to resolve any bubbles 
surrounding the central super-star cluster even though their 
presence in the galaxy is likely. 
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4.3 The Broad Component Gas 

The [O ill] (A5007+A4959)/A4363 intensity ratio can be used 
to determine the electron temperature of the gas from which 
the broad component emission arises. Since, however, [O ill] 
A4363 is only detected as a broad line, this dictates that T e 
estimates using [O ill] can only involve the broad C2 com- 
ponents; moreover these will only be applicable to the in- 
ner core region of the galaxy where A4363 is detected. The 
FWHMs of components C2 in AA4959, 5007 are very similar 
to that of A4363, all being in the range of 320-515 kms -1 
across the inner core region; their surface brightness distri- 
butions are also similar. Also in support of the suitability of 
the [O ill] ratio involving the C2 components is that we do 
not find any indication of C3 or C4 components in the high- 
est S/N emission line profiles of AA4959, 5007 (see Fig G^f) 
for an optimised two-component fit to the A5007 emission 
line from the inner core region). Even so, electron tempera- 
tures derived using only the broad line components are un- 
physically high (> 30,000 K) until densities higher than N e 
— 10 6 cm' 3 are adopted. 

We can attempt to constrain the electron tempera- 
ture and density of the broad line component by using the 
A1663/A4363 and A5007/A4363 intensity ratios, which due 
to their considerably different 2+ upper level excitation 
energies can act as sensitive temperature and density diag- 
nostics between the critical densities of each of these lines. 
The FOS UV spectra presented by TIL96 revealed a weak 
Oiii] A1663 line detection which, due to its high excitation 
energy ought to originate from the same gas emitting the 
broad [Oiii] A4363 that we detect. Using a simulated FOS 
aperture on our broad component flux maps for A5007 and 
A4363, fluxes of 39.6+6.3 x lO" 15 ergs cm' 2 s' 1 for A4363 
and 214.1±14.2xl0' 15 ergs cm" 2 s" 1 for A5007 are pre- 
dicted to fall within the small FOS aperture centred on 
the nucleus of Mrk 996. Together with the A1663 flux of 



■> ; : .- -t 



1.9 



+0.4 



xl0~ 15 ergs cm -2 s _1 measured from the FOS spec- 
trum, these fluxes were then dereddened using an average 
c(H/3) of 0.65 for a simulated FOS aperture on the broad 
component VIMOS c(H/3) map, to yield A4363/A1663 and 
A5007/A4363 intensity ratios of 2.4 and 4.3, respectively. 
The theoretical A4363/A1663 and A5007/A4363 ratios for 
log(iVe) = 5.0-10.0 and T c = 5,000-30,000 K are presented 
in Fig 115b . b and the resultant set of solutions are shown in 
Fig 115b . The two curves in Fig 115b overlap between log(A r e ) 
= 6.5-8.5 and T c = 7,000-14,000 K. We thus adopt T c = 
10,500 ± 3,500 K and log(7V c /cm- 3 ) = 7.25±J;?| for the 
broad component emitting gas. The inferred density range 
is somewhat higher than the value derived from the [Fein] 
A4881/A4658 ratio (Section l4.2[) and supports the conclusion 
that the central regions of Mrk 996 contain very dense zones 
of ionised gas. 



5 CHEMICAL ABUNDANCES 

5.1 Narrow Component Abundances 

Abundance maps relative to H + were created for the N + , 
0+ 2+ , S+, S 2+ , Ar 2 + ions, using the narrow AA6584, 
7320+7330, 5007, 6717+6730, 6312, and 7136-A lines respec- 
tively (via the ABUND task in iraf). Each VIMOS spaxel was 
treated as a distinct 'nebular zone' with its own set of phys- 
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Figure 17. N/H abundance ratio map based on the narrow com- 
ponent emission across the region where [N II] A6584 is seen. It 
should be noted that these abundances are lower limits since an 
upper limit to T c was used for the calculation. 



ical conditions. Ionic argon and sulphur abundance ratios 
were converted to total Ar/H and S/ H using an ICF(Ar) = 
1.87 and ICF(S) = 1.03 respectively (iKingsburgh fc Barlowl 
|l994|). Very small changes were found by adopting the 
ICF prescriptions of Izotov et al. (2006). Since T e maps 
are unattainable for the narrow component gas, we have 
adopted the upper limit of T c — 10,000 K applicable 
throughout Mrk 996 (Section 14. ip . together with the mean 
electron densities derived from the [S n] doublet ratios, of 
170 cm -3 and 10 cm -3 for the inner core and outer re- 
gion, respectively (Table [3}- It is emphasised that since the 
adopted T c for the narrow component gas is an upper limit, 
the abundances derived here are necessarily lower limits. 

In Fig [16] abundance distributions are shown in the di- 
rections from north-east to south-west and south-east to 
north-west that pass through the centre of Mrk 996 (i.e. 
in the X and Y directions across the abundance maps, re- 
spectively, produced by averaging over a three spaxel-wide 
pseudo-slit). The variations in the abundances of 2+ and 
0+ in Fig ma mirror each other in their distribution. As ex- 
pected, the more highly ionised species is dominant in the 
core region, where the UV radiation is expected to be harder. 
It should further be noted that these abundance trends are 
only correct if the temperature is constant throughout the 
gas emitting the narrow line component. 

The O/H elemental ratio was derived by summing 
the + and 2+ ionic abundance maps. Adopting a so- 
lar oxygen abundance of 4.9xl0~ 4 relative to hydrogen 
ijAllende Prieto et al.l 120011 ), we find that Mrk 996 has an 
oxygen abundance of >0.50 solar throughout the inner core 
region and beyond. This is a significantly higher metallicity 
than the ~0.2 solar derived previously by TIL96, which can 
be attributed to the high electron temperature of 15,000 K 
that they employed, based on CLOUDY models. 

An N/H abundance ratio map (Fig I17p was cre- 
ated using an ionisation correction facto r, ICF(N) = 
(0 + +0 2+ )/0+ (|Kingsburgh fc Barlowl Il994h . Again, very 




Figure 15. Broad component T e and N c diagnostic diagrams for (a) 7(A4363)//(A1663), and (b) /(A5007)//(A4363), showing theoretical 
ratios for T c = 5,000—30,000 K and log(7V c /cm — 3 ) = 5.0-10.0. The horizontal solid lines in (a) and (b) denote the observed dereddened 
ratios and their uncertainties (dashed); (c) shows the (T e , JV e ) solutions applicable to Mrk 996 obtained from (a) and (b) and their 
uncertainties. 



s mall changes were fo und by adopting the ICF prescriptions 
of llzotov et al.l (|2006t ). For a constant electron temperature, 
the N/H ratio, shown in Fig 1171 remains constant through- 
out the core at > 0.8 x 10~ 5 but shows a two-fold increase 
north-east of the core region of Mrk 996. A cut across the 
N/O ratio map of Mrk 996, which is much less sensitive to 
the adopted electron temperature than either O/H or N/H, 
shows a similar trend, with log(N/0) — —1.5 in the core, 
rising to ~ —1.2 at a spot north-east of the core (Fig I16[) . 

The mean abundances are listed in Table 5. For the 
revised oxygen metallicity of Mrk 996 the derived N/O, S/O 
and Ar/O ratios fall w ithin the expected range for BCDs 
(e.g. llzotov et al]|2006h . 



5.2 Broad component abundances 

In order to derive nitrogen elemental abundance for the 
broad component gas, we correct for the presence of N 2+ 
using the Nm] A1750 line detected by HST FOS. This high 
excitation energy UV line was not used to correct the narrow 
emission component nitrogen abundances because it ought 



to originate from the high excitation gas responsible for the 
broad line emission. Table [4] lists the ionic abundances de- 
rived for the broad line emission gas derived from optical and 
UV lines. Fluxes are quoted relative to the broad H/3 flux 
within the simulated FOS aperture and are corrected for ex- 
tinction using the average c(H/3) value within that aperture 
in the reddening map derived from the broad components of 
Hq/H/3, H7/H/3 ratios. The physical conditions adopted for 
these calculations are those derived from the broad compo- 
nent emission, i.e. T e = 10,000 ± 3,500 K and log(A" c /cm" 3 ) 
= 7.25 t.Q '75. We consider that the highest allowed T c is 
11,000 K with log^c/cm -3 ) = 7.00 being the lowest allowed 
density, based on Fig 15c. This latter density is consistent 
with the upper density limit inferred from the [Fe ill] diag- 
nostics. For T c = 11,000 K and log(7V c /cm" 3 ) = 7.00 the 
resulting broad-line O/H abundance ratio of 2.3 x 10~ 4 is 
consistent with being the same as for the narrow line emis- 
sion regions (cf. Tables [4] and [5} . As a consistency check, 
we determined the T e for the narrow component VIMOS 
[O ill] A5007 and H/3 fluxes falling within the FOS aper- 
ture that would be required to match the above broad-line 
O/H ratio. This yielded a temperature of ~9,200 K, rein- 
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Figure 16. Spatial variations in narrow component ionic and elemental abundances across Mrk 996 obtained for T e = 10,000 K and N c 
= 170 cm -3 . The available T c is an upper limit so the abundances shown arc lower limits. The patterns should persist if T e is constant 
across the emitting region. Three spaxel-wide averages were made across the IFU aperture in its X and Y directions, corresponding to 
south-east to north-west and north-east to south-west, respectively. 



Table 4. Intensities of UV emission lines from HST FOS spectra and their optical broad component counterparts, measured within a 
FOS aperture superimposed on our broad component flux maps. All intensities are given relative to the broad component F(H/3) = 122.1 
± ll.lxlO - 15 erg s — 1 cm -2 within the simulated aperture. 
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19.2±6.8 
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3.0±1.1 


1.3±0.5 






N/H 


9.4±1.3 


3.8±1.7 


1.7±0.8 


5007 [Oin] 


51.3±6.2 


2 +/H+ 


12.6±1.5 


4.8±0.6 


2.0±0.2 






O/H 


15.7±7.5 


5.5±2.6 


2.3±1.1 



forcing the upper limit of 10,000 K derived in Section |4~T1 for 
the CI emission. S/H abundances were computed using the 
[S ill] A6312 line and an estimated (minor) contribution from 
S + /H + (in parentheses in T able 5) using equation A38 of 
iKingsburgh fc Barlow! fl994). Ar/H abundance ratios were 
derived using the [Arm] A7136 line with ICF(Ar)=1.87 as 
previously. The adopted mean values are listed in Table 5. 

In contrast to the apparent lack of significant variation 
in the abundances of oxygen, sulphur, and argon between 
the two velocity components, it appears that the N/O ratio 
in the broad line region is enriched by a factor of ~20 rel- 
ative to the narrow line region, with the broad line region 



having log(N/0) ~ —0.13. We conclude that the observa- 
tions are consistent with the narrow and broad line regions 
both having an O/H ratio of >0.5 solar. 

5.3 Elemental Helium Abundance 

Narrow and broad component He + /H + abundance maps 
were derived from He I A5876 CI and C2 maps. Abun- 
dances were calcula ted using the Case B He I emissivities of 
IPorter et all (2005) at the respective T c and of the broad 
and narrow emission components. Since these are only valid 
up to log(N c /cm~ 3 ) = 6, a correction factor of 0.9543 was 
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Table 5. Adopted abundances for the narrow (CI) and broad 
(C2) component emission regions throughout Mrk 996. 
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Figure 18. WR 'blue bump' feature summed over the inner core 
of Mrk 996. Overplotted is the combined fi t for 2600 WNL stars 
and 400 WC stars, using the templates of ICrow thcr & H adfieldl 
l|2006h . 



applied to the broad component emissivity to extrapolate 
up to log(-/V e /cm _3 ) = 7.25. At these densities, collisional 
excitation from the He 2s 3 S metastable level by electron 
impacts contributes to the observed br oad He I A58 7 6 com - 
ponent; this effect is accounted for bv lPorter et al.l ((2005) . 
Fig 1191 shows X and Y direction cuts across the broad and 
narrow component elemental He/H maps. We have assumed 
that the He 2+ /H + fraction is negligible. Good agreement be- 
tween the spatial variations in the He/H abundance for each 
component is observed. No relative He/H enrichment is seen 
across the components; within the uncertain ties both fall be- 
tween the pre-galactic He/H ratio of 0.08 (fLuridiana et al.l 
120031 ) and the solar ratio 0.10 (|Loddersll2003l ). This contrasts 
with the localized N/H enrichment north-east from the nu- 
cleus that spatially correlates with a peak in EW(Rf3) (see 
Section [6]). The slight decrease in He/H in the core region 
might be attributable to the presence of He 2+ . We find this 
unlikely, however, since no nebular He II A4686 emission has 
been detected (non e was reported in the recent study by 
iThuan et al.l ([2008)), and in addition, other similarly high 
excitation lines such as [Ariv] are not present in the spec- 
trum of the galaxy. The derived He/H abundances at the 
centre of Mrk 996 agree well with those of TIL96 who cal- 
culated He+/H+~0.09 from the He I A5876 emission within 
their FOS aperture. 

Table [5] summarises the derived abundance for the nar- 
row and broad line emission regions of the galaxy. 
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a Upper limit adopted T so that CI abun- 
dances are lower limits. 

b Applicable to a central 0."86 (FOS) aper- 
ture; see the text for details. 
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Figure 19. Spatial variations in the narrow (blue and red) 
and broad (green) component He/H elemental abundance across 
Mrk 996. Three spaxel-wide averages were made across the IFU 
aperture in its X and Y directions, corresponding to the south- 
east to north-west and north-east to south-west, respectively. 



6 WOLF-RAYET STARS AND STELLAR AGE 

Wolf-Rayet (WR) features were identified in Mrk 996 by 
TIL96 in their single aperture FOS spectroscopy of a 0."86 
region at the centre of Mrk 996. The VIMOS IFU spec- 
tra also show a broad Wolf-Rayet (WR) stellar feature at 
4650 A, attributable to a mixture of late-type WN (WNL) 
and WC stars at He II A4686 and N ill A4640 (with an ad- 
ditional weak WC stellar C IV feature at 5808 A). The WR 
spectral signatures are seen to extend throughout the core 
region of Mrk 996 and show excellent agreement with the 
spatial extent of the broad [O III] A4363 emission (Fig [8}. A 
map of the 4650 A WR feature is shown in Fig 1201 this has 



been decontaminated from the emission by nebular [Fe III] 
A4658 by fitting a Gaussian profile to the line and subtract- 
ing its flux from the integrated flux of the 4650 A feature. 
After summing the spect ra over the core regio n , LM C WR 
spectral templates from ICrowther fe Hadfieldl ([2006) were 
used to fit both spectral features. Although the WC feature 
is weak, its contribution is essential in fitting the broad wings 
of the 4650 A feature. At the distance of Mrk 996 (22 Mpc) 
the two WR features indicate the presence of ~2600 WNL 
stars and ~400 WC stars. Assuming a similar depth for the 
area over which these spectral features are summed, 1.7x2.3 
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Figure 20. Emission map of the blue WR feature (created by 
integrating over the full emission feature and removing contam- 
inant flux from [Fein] A4658), showing the distribution of WR 
stars in the inner core region of Mrk 996. Overlaid are the Ha 
narrow component flux contours from Fig [l] 



arcsec 2 (4.6 xlO 4 pc 2 



results in a volume density of 3 x 10 -4 
WR stars per pc a in the core region. 

We can compare our WR number estimates with those 
of TIL96, who estimated that 600 late- type WN stars and 74 
WC stars were responsible for the WR emission seen within 
their 0."86 diameter FOS aperture. Scaling our own core 
region WR number estimates for the difference in aperture 
size would predict 390 WNL stars and 60 WC stars within 
the FOS aperture. While the agreement between the two 
estimates is quite good, the decrease in WNL star numbers 
for our scaled estimate can be attributed to the non uniform 
density of WR stars in the core region (see Fig I20[) . In con- 
trast, the WC 5808 A feature is very centrally peaked, thus 
giving a good agreement between the two estimates. 

Luminosities of hydrogen recombination lines, particu- 
larly H/3, can provide estimates of the io nizing flux present, 
assum ing a radiation-bounded nebula (jSchaerer fc Vaccal 
1 19981 ). Thus, the equivalent width (EW) of H/3 is commonly 
used as a stellar age indicator at a given metallicity. A map 
of EW(R(3) across Mrk 996 (derived from integrating over 
the full emission line profile) is shown in Fig [2T]and displays 
a surprising morphology compared to the central location of 
the ionizing sources in Mrk 996. A peak is seen directly NE of 
the centre of the galaxy, which correlates well with the peak 
in [Nil] A6584 emission shown in Fig 1171 We can use the 
EWQ10) map, in conjunction with the metallicity map (de- 
scribed in Section [5. ip . to estimate the distribution of stel- 
lar ages throughout Mrk 996. We find that the core region 
EWQ10) corresponds to a stellar age of 4.5 Myr whereas 
the £W(H/3) peak NE of the core corres ponds to a slightly 
youn ger stellar age of ~3 Myr (cf. fig 7 of ISchaerer fc Vaccal 
I1998T ); these ages are sufficiently long to allow for the pres- 
ence of WR stars in the galaxy. 

An estimation of the number of O-type stars was made, 
assuming that all ionizing photons, Qo bs 5 are produced by 



O and WR stars and utilising the H/3 luminosity, L(H/3), 
integrated over the whole galaxy (i.e. using the dereddened 
combined CI and C2 H/3 flux from the inner and outer re- 
gions, and at a distance of 22.3 M pc). The absolute numb er 
of O stars was derived using (e.g. lFernandes et all l|2004h ) 
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where QJj bs is the emission rate of ionising photons, as 
derived from L(H/3), and Qwr, Qo7v are the emission 
rates of ionizing photons by WR stars and 07V stars, 
respectively, with TVwr being the total number of WR 
stars. rio(t) is a conversion parameter for the propor- 
tion of 07V stars relative to all OV stars a nd is es- 
timated as being >0.2 by ISchaerer fc Vaccal (|l998h at 
the metallicity of Mrk 996 (~0.5Z©). Adopting rio(t) 
=0.2, Q W r = Qorv = 1.0 x 10 49 s _1 l|Schaerer et all 
I1999T ). and Qg bs (Cl+C2)=3.36x 10 53 photons s" 1 we find 
iVo < 153,000 yielding N WR /N v > 0.02. As a lower 
limit, this ratio is in good agreement with those pre- 
dicted by evolutionary synthesis models by previous stud- 
ies; at the predicted age of < 5 Myr for the Mrk 996 
starburst, ICervifio fc Mas-Hessd (1 19941 ) find that this 
Nwr/Nqv fractio n is typ ical at metallicities of 0.4 Zq, while 
ISchaerer fc Vaccal l| 19981 ) predict N W r/(N wr + No) = 0.04 
for a metallicity of 0.25 Zq. 

Using the H q-SFR (star formation rate) relation of 
iKennicuttj l| 19981 ) and the total Ha luminosity (C1+C2) 
of 4.83 xlO 41 erg s _1 we obtain a global current SFR of 
3.8 Mq yr -1 . Taking into account the subsolar metallicity 
of M rk 996 this can be reduced to 2.7 M yr" 1 (|Lee et all 
2002). However, in the event that the C2 component emis- 
sion is partly attributable to shock-excited gas (see Section 
I7.1[l so that it cannot be included in the estimate of Qo bs , a 
SFR lower limit of 1.3 Mq yr -1 is derived based on Qo(Cl) 
= 1.31 xlO 52 photons s" 1 alone. 



7 DISCUSSION 

7.1 Diagnosing the ionisation mechanisms 

In order to gain insight into the mechanisms responsible 
for the multi-component line emission fro m Mrk 996, we 
made use of the classic diagnostic diagrams of lBaldwin et al.l 
(1981) (the BPT diagrams). These are employed to clas- 
sify galaxies according to the dominant excitation mecha- 
nism of their emission lines, i.e. either photoionisation by 
massive stars within H II regions or photoionisation by 
non-thermal continua from active galactic nuclei (AGN). 
The diagrams consist of excitation-dependent, extinction- 
independent line ratios: log([Olll] A5007/H/3) versus ei- 
ther log([Nll] A6584/Ha) or log([S ll](A6716,31)/Ha). Star- 
forming galaxies fall into the lower left region of the diagram, 
AGN host galaxies fall into the upper right region and Low- 
ionisation Emission Line Regions (LINERs) fall in the lower 
right region. The separation is not as clear, howev er, for 
low metallicity AGNs (see e.g. IStasiriska et al.l ([2006)), only 
a handful of which ha ve been proposed to exist thus far 
jlzotov fc Thuanll2008l ). 
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Figure 21. Map of the equivalent width of H/3 across Mrk 996, 
derived from the full emission line profile. Contours range from 
70-180 A in steps of 30 A. 



iKewlev et" al. (2001) calculated the first starburst grids 
that attempted to match optical diagnostic diagrams based 
on purely empirical data, by coupling the mappings III pho- 
toionisation code and STARBURST99 population synthesis 
models. Whereas previous semi-empirical studies involved 
only solar-metallicity calculations, their line ratios were 
computed for a range of metallicities (Z = 0.05-3.0 Zq) and 
ionisation parameters (q = 5 x 10 6 -3 x 10 s cm s _1 ), where q 
is the maximum velocity of an ionisation front that can be 
driven by the local radiation field and is relate d to the non- 
dime nsional ionisation par ameter U = q/c, iDopita et al.l 
120011 ). Utilizing these grids, IKewlev et all i|200l] ') determined 
a 'maximum starburst line' (shown in Fig I22[) , above which 
the flux ratios of an object cannot be fitted by pure star- 
burst models alone. Ratios lying above this photoionisa- 
tion boundary require additional sources of excitation such 
as shocks or AGNs. As mentioned above, this picture has 
recently been revised and it is now believed that law-Z 
AGNs, should they exist, would occupy similar regions on 
the BPT diagram as normal H II galaxies; based on an anal- 
ysis of a large sample of Sloan Digital Sky Survey galaxies, 
IStasihska et al.l (|2006T ) have found that composite AGN/H n 
region models with Z < 0.4 Z© he very close to and below 
the pure H II region sequence, whereas when considering also 
the regi me of Z > 0.6 Zp th ey established that the locus be- 
low the IKewlev et all (|200ll ) line allows for an AGN-excited 
gas component of up to 20 per cent. 

It is therefore instructi ve to see wha t area of the BPT 
diagram Mrk 996 occupies. IKewlev et al.l l|200lf ) found that 
since the [S n] (6716,31) /Ha ratio is affected by the density of 
the line emitting region and it is therefore less accurate as a 
diagnostic. Hence considering the large density difference be- 
tween the narrow and broad component gas within Mrk 996 



we only considered the [O ill] A5007/H/3 vs. [N n] A6584/Ha 
ratios (Fig [22} . Since no broad component is detected for 
[Nil] A6584 we have not been able to create separate broad 
and narrow diagnostic diagrams and instead present (i) a 
narrow component diagram (Fig [22] top panel), and (ii) ra- 
tios derived from the entire emission line profiles (i.e. inte- 
grated over broad and narrow component emission, where 
both exist; Fig [22] bottom panel). Spaxels corresponding to 
the core region of the galaxy are shown as red crosses, with 
outer region spaxels plotted as green crosses. Also shown is 
the maximum starburst line of IKewlev et al.l l|200lf) . 

The data points for the outer narrow line region of 
Mrk 996 mostly straddle the 'theore tical' upper limit for 
pure photoionisation. Based on the IKewlev et al.l (|200ll ) 
grids, the flux ratios lie within a metallicity range of Z — 
0.3-0.5 and along and above an ionisation parameter track 
for q — 3 x10 s (the maximum ionisation parameter com- 
puted). The narrow component emission in the core region 
occupies the top left of the distribution, which corresponds 
to a higher ionisation parameter and lower metallicity in the 
starburst grids, while a significant number of data points 
fall outside the SFG (star-forming galaxies) area of the plot. 
When the broad component is included in the diagnostics, 
the location of the outer region spaxels move to the left, 
whereas the core region spaxels move downwards and to the 
left for both ratios to occupy a lower ionisation parame- 
ter and metallicity region. This is most likely due to the 
suppression of the [O in] A5007 broad component emission 
within the high density core region, which would lower the 
combined narrow plus broad [O ill] A5007/H/3 ratio. 

As previously mentioned, exceeding the 'maximum star- 
burst' line can also be indicative of a contribution to the 
ionisation by fast shocks. However, the amount by which 
this line can be crossed before shocks must be a substantial 
source of ionisation has not b een established thus far. Sim- 
ilar diagrams were created bv lCalzetti et all (|2004T) for four 
starburst (SB) galaxies using ratios from HST WFPC2 im- 
ages, such as [O in] A5007/H/3 vs. [S n] A6716,31/Ha. They 
found that a number of regions in the SB galaxies lay above 
the maximum starburst line and as a result defined regions 
of non-stellar excitation as either lying above and to the 
right of it or where the [Sll]/Ha ratio is high enough to 
be compatible with non-stellar sources, i.e. [S ill /Ha > 0.7- 
0.9 l|Veilleux fc Osterbrockl 19871 ; IShull fe McKeej|l979T ). The 
log([Sll] A6716,31/Ha) ratio is traditionally used as a diag- 
nostic for shock-excited gas because shock models predict 
that relatively cool high-density regions shou ld form behin d 
the shock front and emit strongly in [Sn] (|Dopital Il978h . 
The [S ll]/Ha flux ratios for the narrow component emission 
across Mrk 996 lie between 0.03-6 x 10 -4 , well below the high 
[S ii] /Ha rat i os pre dicted for shocks. On the other hand, 
iThuan et all |2008i) employed a model with shocks at a ve- 
locity of 250 km s _1 that reproduced the observed intensity 
of the high-ionisation potential [O iv] 25.9 fim line detected 
in the Spitzer spectrum of Mrk 996. In conclusion, Fig I22l bv 
itself cannot provide firm evidence for or against the pres- 
ence of AGN activity within Mrk 996 given the galaxy's 
rather low metallicity status, even though the integrated 
broad component Ha luminosity of 2.5xl0 41 erg s" 1 this at 
the lower limit of those measured from rare low-Z BCDs 
suspected of h arbouring active galac tic nuclei (3xl0 41 - 
2xl0 42 erg s" 1 llzotov fc Thuanll200etl . Although photoion- 
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isation appears to be the dominant excitation mechanism, 
we cannot rule out a contribution from shocks, particularly 
with respect to the broad component emission. We draw at- 
tention to the significant bias inherent in BPT diagrams of 
galaxies whose analysis does not involve a separate consid- 
eration of narrow/broad line diagnostics. 

7.2 What is the origin of the broad line emission? 

Mrk 996 is not the only BCD t hat shows evidence for large 
line widths. Ilzotov et alj l|2007h presented an analysis of a 
large number of low metallicity BCDs exhibiting broad emis- 
sion lines with inferred gas expansion velocities of 1000-7000 
km/s. WR stellar winds have been proposed as a possible 
mechanism for producing such broad line widths, but thus 
far no correlation between the presence of b road nebular 
lines and broad WR features have been found (Ilzotov et al.l 
120071 ). We find that in Mrk 996 the surface brightness of 
the 4650A WR feature does not correlate tightly with that 
of [O in] A4363 (Fig 23), further indicating that the mech- 
anism responsible for the broad line emission is not active 
onl y in the reg i ons where W R stars are present. 

iRov et aH (|l99ll . Il992h investigated a variety of mech- 
anisms to explain broad nebular gas components in H II 
regions, including electron scattering, stellar winds, super- 
nova remnants and superbubble blowouts. However, each 
of these mechanisms was deemed unsatisfactory. An explo- 
ration of blowout mecha nisms within low-me t allicit y H II 
regions was conducted by iTenorio-Tagle et al.l |l997h using 
a hydrodynamical calculations to try to match a sample of 
low-metallicity H II regions with broad emission line compo- 
nents. They proposed that the absence of strong radiative 
cooling in the low-metallicity ISM within these regions could 
delay the action of Rayleigh- Taylor (RT) instabilities that 
can fragment the expanding shell within a blowout. How- 
ever, this mechanism may not be applicable to Mrk 996, as 
one of the required conditions is a low-density environment, 
whereas we derive high densities for its broad component 
gas (log(JV c )~7). 

Broad line emission (FWHM w 400 km/s; C2) is de- 
tected throughout the central regions of Mrk 996, through- 
out the bright starburst region and out to >5". Close to 
the nucleus of the galaxy the broad Balmer line compo- 
nent splits into two additional broad components (FWHM 



200-300 km s" 



which trace the kinematics of a two- 



armed mini spiral, but component C2 does not show evi- 
dence of velocity splitting at our level of spectral resolution 
farther out (Figs 3 and 4). Given the youth of the star- 
burst region in Mrk 996 (<5 Myr) containing large numbers 
of young O and WR-type stars (Section 6), and presence 
of very dense zones of gas (Section 4.3), it is possible that 
the broad line emission partly originates from a turbulent 
mixing layer forming on the surface of dense cool clumps 
which are subjected to irradiation and hydrodynamic abla- 
tion from the hot wi n ds from young mas si ve sta r clusters 
l|Begelman fc Fabianl (|l990h : ISlavin et al] (|l993l )). Recent 
work on the nearby NGC 1569 and M82 has revealed that 
the broad line component observed in those starburst galax- 
ies can be explained as originating wi thin turbulent layers 
on the surface of the dense gas clumps ( Wcstmoquctto et al.l 
l2007al lbTl. Motivated by this hypothesis. lBinette et all (|2009l 'l 
have successfully modelled the broad underlying components 



seen in the Balmer and [O ill] lines of NGC 2363 via the 
inclusion of a turbulent mixing layer and conclude that the 
broad profile results from radial acceleration of photoionised 
turbulent gas. Their models also correctly predict the ab- 
sence of a broad component in [S n] and [N n] lines. Al- 
though the applicability of these results to Mrk 996 may be 
limited due to the high densities of Mrk 996 's broad compo- 
nents (JVc =100 cm" 3 is adopted for NGC 2363) and their 
much broader predicted line widths (±3500 kms" 1 ), their 
results lend a great deal of weight to the turbulent mixing 
layer hypothesis. 

This mechanism could also be at work here. It is prob- 
able that the central region of Mrk 996 consists of many 
small, dense clouds of relatively low volume filling factor 
dispersed within the young star clusters of the starburst. 
This results in a large cloud surface area with which the 
copious ionizing photons and fast stellar winds can inter- 
act. If the presence of broad emission indicates that strong 
wind-clump interactions are taking place, then by extension, 
material from these interaction sites must be being stripped 
off and entrained into the cluster wind flows contributing to 
the overall appearance of component C2. 



8 SUMMARY AND CONCLUSIONS 

Mrk 996 is a complex system containing an extended star- 
burst region (216-432 pc in radius) with a centrally con- 
centrated Wolf-Rayet population in the nuclear regions 
(<200 pc). The young stellar population in the nuclear 
super-cluster contains ~3000 WR stars and < 150,000 07V- 
type equivalent stars. The age of the nuclear starburst is 
estimated to be ~4.5 Myr but there is evidence that the 
extended starburst is younger by ~1 Myr. The presence of 
an old stellar population (>1 Gyr) has been independently 
established (Thuan, Hunt & Izotov 2008). The current star 
formation rate is ~2Mg yr -1 . The spatially mapped ion- 
ized gas shows a composite emission spectrum consisting of 
narrow and broad lines. The Balmer line velocity structure 
close to the nucleus confirms the presence of a spiral struc- 
ture confined in the inner 160 pc. The Ha broad compo- 
nent has a large integrated luminosity of 2.5xl0 41 erg s _1 ; 
this is at the lower limit of those measured from the rare 
low- Z BCDs suspected of harbouring active galactic nu- 
clei (|lzotov fc Thuanll2008h . Our VIMOS IFU analysis has 
enabled a separate analysis of the physical conditions and 
chemical composition of the narrow and broad emission line 
regions. The broad line region in the nuclear starburst is 
very dense (10 7 cm -3 ) whereas the narrow line component 
is of lower density (<10 3 cm -3 ). 

The upwards revised oxygen metallicity of Mrk 996 is 
> 0.5 Z Q (12 + log O/H w 8.37). The broad line region 
is nitrogen-enriched compared to the narrow line region by 
<1.3 dex. However, no relative enrichment of He, O, S, and 
Ar is inferred. The S/O and Ar/O abundance ratios in the 
two components, as well as the N/O ratio in the narrow 
line component, are typical of those in H H galaxies and 
dwar f irregulars (e.g. Ivan Zee fc Havnesll2006l ; Ilzotov et al.l 
2006). The narrow- line N/O ratio in particular is exactly the 
value expected based on the galaxy's colour, and follows the 
metallicity- lumin osity relationship for isol ated dwarf galax- 
ies (cf. figs 8, 9 of Ivan Zee fc Havnej200^ for M B = -16.7, 
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Figure 22. Emission line diagnostic diagram for Mrk 996. Each data point represents a spaxel in the ratioed dereddened flux ma ps 
corresponding to [Om] A5007/H/3 and [Nil] A6584/Ha. The empirical diagnostic 'maximum starburst' line from lKewlev et"aL is 
also shown, as are the positions of emission line ratios for star- forming galaxies (SFG) and active galactic nuclei (AGN). The top panel 
represents narrow component emission line map ratios whereas the bottom panel shows ratios derived from the full emission line profile 
(i.e. narrow + broad component emission, where detected). 



B — V — 0.44, TIL96); its mean current value can be inter- 
preted as the result of the slow release of nitrogen from the 
intermediate-mass stellar population over the last few Gyrs. 
On the other hand, the high N/O ratio in the broad line 
region of the inner galaxy is consistent with the presence of 
numerous evolved massive stars (e.g. WNL-type and Lumi- 
nous Blue Variables) and can be attributed to the cumula- 
tive effect of their N-enriched winds. An elevated N/O could 
not in this case be due to the outflow of oxygen-enriched gas 
from supernova explosions: the absence of elevated (or sup- 
pressed) O/H, S/H and Ar/H ratios in the broad line gas 
implies that recent supernova ejecta are probably not im- 
plicated in its excitation; at the young age of the Mrk 996 
starburst only a few very massive stars (>5OM0) would 



have exploded as supernovae (e.g. IWooslev et al ] |2002h . The 
mild two-fold increase over the mean of the narrow-line N /H 
(Fig I19[) spatially correlates with a local peak in EWQ10) 
in a region (Fig I2ip where no WR stars are seen. This could 
mean that even though the narrow-line gas is fairly well 
mixed, some localized N-enrichment has already occurred in 
an area dominated by normal OB-type stars which is slightly 
younger than the nuclear starburst. 

Finally, we would like to draw attention to the nature 
of the broad [O ill] A4363 line detected in Mrk 996. Even 
though this galaxy represents an extreme case in harbour- 
ing a very dense ionized component which dominates the 
excitation of this line, it is by no means unique (see e.g. 
Ilzotov et ai1l2007l ). In such cases where a substantial frac- 
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Figure 23. Grey scale and contour emission line maps of the Hell A4686 WR feature (left panel) and nebular [O in] 4363 (right panel). 
Their contours suggest that the ionisation mechanism responsible for the broad emission line components (represented by [O III] A4363) 
may not be strictly due to the WR stars. Flux contours are in units of xl0 _15 erg s _1 cm -2 per 0.33x0.33 arcsec spaxel. 



tion of the [O ill] 4363 and 4959, 5007 A lines arise in differ- 
ent gas components, the integrated auroral to nebular ratio 
would not be representative of the electron temperature of 
the overall H n region; using it would result in biased physi- 
cal conditions and chemical abundances. The problem would 
be aggravated in analyses based on low dispersion, low spa- 
tial resolution spectra and may therefore have implications 
for t he so-called 'heating problem' of blue compa ct galaxies 
(e.g. [stasiriska fe Schaererlll999l ; IPequignotJkoOSl ). 
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